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Abstract 

t— i ; 

' The upper bound on the mass of the lightest Higgs boson is provided in 

the supersymmetric standard model with multi-Higgs doublets, up to two- 
ff) ' loop order. Relatively large corrections are expected from the experimentally 

^ , unconstrained extra Yukawa couplings. We calculate it including the two-loop 

C****" ■ leading-log contributions as a function of the scale A, below which the theory 

remains perturbative. Even for A = 10 4 GeV, we obtain the upper bound of 
! 140 GeV, which is heavier than that of the MSSM by 10 GeV. 

m 
O . 

1 Introduction 

\ There is a famous upper bound on the mass of the lightest Higgs boson in the 

minimal supersymmetric standard model (MSSM), i.e., which is lighter than Z 
boson at tree level |l| . This comes from the fact that the quartic couplings of the 
neutral Higgses are determined by only the weak gauge couplings g and g' due to the 
supersymmetry (SUSY). However, since the top Yukawa coupling is large, this tree- 
level upper bound receives large corrections ^ ||, 0, p[. Further, QCD correction, 
which begins from a two-loop order, is also large M 0, |8|, [| [l0| [TTJ , because the 
$_J ■ strong gauge coupling a s first appears at two-loop order, in other words, two-loop 

is the lowest order in terms of a s . Consequently, the upper bound can be set about 
130 GeV in the MSSM when all the squark masses are at order of 1 TeV. 

Here we will provide the upper bound on the mass of the lightest Higgs boson 
in the supersymmetric standard model with arbitrary number of Higgs doublets, 
which is the minimal extension of the MSSM in the meaning that its Higgs sector 
is composed of doublets only. This model has the same upper bound as the MSSM 
up to one-loop level. So, we must include at least two-loop corrections to the 
upper bound in order to see the difference from the MSSM. Actually, the two-loop 
order corrections include the contributions of the extra Yukawa couplings, which 
are characteristic of the multi-Higgs-doublet model and can be relatively large. 

We demand that the theory remains perturbative up to the scale A in order 
to set the upper bound on the extra Yukawa couplings, which are unconstrained 
experimentally. Then the resulting upper bound mh depends on the scale A, and 
becomes more stringent as A is larger. For example, we obtain nih — 128 GeV in 
the case of A = 10 16 GeV, which is almost the same as the upper bound in the 
MSSM. On the other hand, for A = 10 4 GeV, which is one order above the SUSY- 
breaking scale ms, the upper bound is lifted up by 10 GeV due to the large extra 
Yukawa couplings. Throughout this letter, the soft SUSY-breaking squark masses 
are assumed to be 1 TeV and degenerate for simplicity. 
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2 Up to one-loop order 



We denote the number of the Higgs doublets as 2n. |] Then the superpotential of 
this model is written as 



W = ^2 h2iQ 3 H 2 iU 3c — ^2 h2i-iQ 3 H 2 i- 1 D 3 



■)3c 

1=1 i=l 



where 



+ X! / J, 2i-i,2jH 2 i-iH 2 j + h.c. H , (1) 

i J = l 

where Q 3 , C/ 3c and D 3c are the superfields of the left-handed quark doublet in the 
third generation, right-handed top quark and right-handed bottom quark respec- 
tively, and H 2 i— i and H 2 i are the superfields of the Higgs doublets. The ellipsis 
denotes the terms involving quark superfields in the first and second generations 
and lepton superfields, which are irrelevant to our result. 

With general soft SUSY- breaking terms, the tree-level Higgs potential of this 
model is given by 

+4>°JmI^ + <$Ml<j>l - (VSmJ^S + h.c), 

(2) 

J , 0° (3) 

are the neutral components of the Higgs doublets, which couple to the down-type 
quarks and to the up-type quarks respectively. The matrices M 2 and M 2 are n x n 
Hermitian matrices and M 2 U is generally annxn complex matrices. The couplings g 
and g' are the gauge couplings of SU(2)l and U(1)y respectively. We concentrate 
our attention only on the neutral Higgs fields since only neutral fields may take 
VEVs after the breakdown of the electroweak symmetry. 
Now we can rotate the Higgs fields so that 

«) = = (4') = ■ = (4n) = 0, (v d ,v u : real), (4) 

where prime denotes a quantity after the field transformation. 

Just like in the MSSM case, we regard a diagonal element of the mass squared 
matrix, (</>|M 2 |</>), as the desired upper bound. Here the matrix M 2 is the mass 
squared matrix of the neutral Higgses, and the field <fi is the real part of the field <fi° 
defined by 

0° W cos/3 sin/3 , taaB= Vu ) (5) 

X° ) -{ -sm/3 cos/3 ) { <j% J' ^ ~ vj (5) 

and 

<tP = ±={<t> + vp). (6) 
Then we can restrict the effective potential to tp®' and 4>2 directions. 

+mf\^\ 2 + m?|$T - K 2 2 <^' + h.c). (7) 



1 The number of the Higgs doublets must be even because of the anomaly cancellation. 
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This is the same form as the Higgs potential of the MSSM. Therefore it is a trivial 
result that the desired tree-level upper bound is the same one as that of the 
MSSM, i.e., m h = m z cos/3 @. 

Next we consider the one-loop radiative corrections to this tree-level upper 
bound. The one-loop corrections to the effective potential mainly come from the 
loops of quarks and squarks in the third generation because of the large sizes of 
their Yukawa couplings. In the MS scheme, this is given by 



32tt 2 




^ qV Q 2 2' ^ qy Q 2 2 

q=hte,biM " 1=t,b 

(8) 

where rriq and m q are Higgs-ficld-dependent masses of the squarks and quarks re- 
spectively, and defined by 

m 2 + m 2 ±|%A t </>° +%/i^o*| 5 (9) 

--m 2 + ml±\ t h b A b (/> d + ^n J 'h b \, (10) 

m 2 = \%<I> U \ 2 , (11) 
m 2 b = Wf, (12) 

m 2 z = \{g 2 +9 l2 ){\cp Q d \ 2 + \^ u \ 2 ) , tan/3 = |^|/|^|. (13) 

Here h t and h b are the row vectors of the Yukawa couplings defined by 

h t ee t (h 2 ,h 4 , ■ ■ ■ ,h 2n ) , h b ee t (h 1 ,h 3) ■ ■ ■ ,h 2n -i) , (14) 

and A t = diag(A 2 , A 4 , • • • , A 2n ), A b = diag(Ai, A s , • • • , A 2n _i) are the soft SUSY- 
breaking A-parameter matrices, and /i is the /i-parameter matrix whose 
component is fJt2i—i,2j- Here we assume the soft SUSY-breaking squark masses 
to be degenerate and take the common mass m. We also neglect the small D-term 
contributions for simplicity, but including these contributions does not change the 
following discussion. 

We should calculate the second derivative of = + V^ lloop ^ in terms of 
(j> to obtain the upper bound nih- In the basis (||), the effective potential restricted 
to the directions of cf)®' and 4» 2 again becomes the same form as that of the MSSM. 
Hence we can easily see that the upper bound is also the same as the one in the 
MSSM case and conclude that m,h = 140 ~ 150 GeV, which is realized in the case 
of the large left-right stop mixing. 

Thus m,h is the same as the one in the MSSM up to one-loop order. At two-loop 
order, however, the extra Yukawa couplings first contribute to m/, as well as the 
strong gauge coupling a s . Then since two- loop order is the lowest order in terms of 
these couplings and one-loop corrections are the same order as the tree-level value, 
the two-loop corrections might be relatively large. In fact, the QCD correction 
lowers the upper bound by more than 10 GeV in the MSSM due to the large size of 
a s . Since the extra Yukawa couplings are unconstrained experimentally, we will set 
their upper bounds by requiring that the theory remains perturbative up to some 
scale A. The maximal values of the extra Yukawa couplings under this requirement 
might be large, so that the two-loop corrections to due to these couplings are 
as large as the QCD correction. So in this letter, the upper bound mh is calculated 
using the RGE improved effective potential approach in the two-loop leading-log 
approximation in order to see how different the upper bound mh is from that of the 
MSSM. 
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3 Upper bound formula 



For simplicity, we assume that the squarks are degenerate and their masses are 
1 TeV throughout this letter and neglect 0(g 4 ,g 2 g' 2 ,g' A ) contributions. We can 
obtain the upper bound formula in a quite similar way to Ref.|l3|, in which the 
analytic expression of the lightest Higgs mass is provided in two-loop leading-log 
approximation in the MSSM. 
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where /3 mt is defined as 



Q 



9 (/i' 2 sin /3) _ h' 2 sin /3 



16?r 2 



(15) 



(16) 



and (3 mt [m f ] stands for /3 mt defined at the scale Q = Wit- The renormalization scale 
is set to be nit, and rh~z = mz(m t ), rn~b = mb(m t ). The mass parameter m t is the 
on-shell running mass m t (M t ), where M t is the pole mass, and the relation between 
rh~t and M t in the MS scheme is 

— - M% 

mt ~ l + 4a s (M t )/37r' ^ ' 

The parameter X t is the mixing parameter of the left-right stop mixing, which is 
defined by X t = \A' 2 + // 12 cot/3|, and mg = y/ m 2 + m 2 can be regarded as the 
SUSY-breaking scale. 

The factor m 2 z cos 2 2(3 in the first line of (jlj) corresponds to the tree-level contri- 
bution and the remaining factor of the first line represents the anomalous dimension 
of the lightest Higgs field between the SUSY-breaking scale ms and the renormaliza- 
tion scale mt. The second line of the formula is the part of the one- loop correction, 
and the rest is the two-loop leading-log contributions. 

It is clear from the above expression that the upper bound m^ gives the absolute 
upper bound when tan/3 is large and (3 mt \m t ] has its maximal value. The (3- function 
(3 mt is expressed by using the step functions as 

(3 mt = /3v +cos 2 /3( 7 «- 7 W) 



75 s 2 (6 



'GQ 



'Gil) + ^2 



\h' 2 f 



2 ^ 

i=l 



1 ™ 



29 7, r\ 



'4>'u 



(18) 



i=l 



where g s is the strong gauge coupling and the ellipsis denotes the terms involving 
g, g' ', • • •, which can be neglected because they are irrelevant within our approxima- 
tions. A step function Op is defined as one above the mass of the particle P and 
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zero below it, and Bab = @a9b- (G denotes gluino, Q, U and D are the squarks in 
the third generation, <(/• are Higgses and <^ are Higgsinos.) Now since we consider 
below the scale Q = m{= uiq = rrifj = rag), the maximal value of j3 mt is 

o n , n 

0m t =-8 5s 2 + 3/i; 2 sin 2 /3+-^|/ l ' 2l | 2 + -^|^. 1 | 2 , (19) 

i=l i=l 

which is realized in the limit that all the Higgs bosons are lighter than the top 
quark. 

Here we rotate in the Higgs field space again and make the only one pair of the 
Higgs doublets couple to quarks. The top and bottom Yukawa couplings in this 
basis are 

n n 

h? = ^2\h' 2i \ 2 , c = £i>4-ii 2 , (20) 

i=l i=l 

respectively. (The quantities in this basis are denoted by the double primed ones.) 

However, this limit does not coincide with the limit that the lightest mass eigen- 
value saturates the upper bound m^. In fact, not all the Higgses have to be heavy 
to saturate m^, but some of them may be relatively light unless they mix with the 
lightest mass eigenstate. In this case the scale at which (3 mt is estimated is not 
above all of the Higgs mass eigenvalues and thus (3 mt [mt] that should be used in 
the upper bound formula ( |l5| ) is smaller than the maximal value (|l^ ) by the con- 
tributions of the "extra" Yukawa couplings corresponding to the "extra" Higgses 
that are decoupled. Here we denote the word "extra" as the meaning of giving no 
masses to the quarks. For example, the extra Yukawa couplings mean h! 2i and h' 2i _ 1 
(i > 2). The absolute upper bound is thus realized when all Higgses that do not 
mix with the lightest one are light and the extra Yukawa couplings that correspond 
to heavy extra Higgses are small enough so that (3 mt is close to ( |l9| ) if the fixed 
values of h' 2 ' and h'{ are given. Thus, from now on, we will restrict the parameter 
space to this region. 

Unlike the MSSM case, there is no upper bound on h 2 and h'{ due to the 
existence of unconstrained coupling constants h! 2i and h' 2i _ 1 (i > 2). Then we put 
the requirement : 

"The theory remains perturbative up to some scale A." 

By substituting ( |l9| ) for /3 mt [mt] in the formula (fi"5|)and setting tan/3 large, the 
desired upper bound is obtained as 

3mf ln m s 



2tt 2 v 2 mt 
mi f 1 Xf M Xl . , m| Xf , , m s 



- ' ' 2 m s 12m| 8m 2 m s m t 



1 

"8^ 



.3 » 2 1 » 2 ,X 2 X 2 . m s 

( + ^i max - 32^)^(1 - ^) ln — 

gg| x l ln m s 

2v 2 m s m t 



m t , 3 , "2 , 1,"2 oo \n ms \ 2 

— + o"2max + ^"lrnax ~ 327Ta s )(ln — 

v £ 2 I m t 



(21) 



where h 2max and h'{ majx are the values of h 2 and h![ that maximize the combination 

3," " '" 



\h 2 2 + 5^i 2 under the above requirement. 
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4 Maximal values of h 2 ' and h'{ 

Now we will compute ft^'max an d h" max . First, we should notice that the RGEs of 
h'2 and h" are the same as that of h t and hb in the MSSM because h' 2 ' and h" are 
the sole pair of the Yukawa couplings that couple to the up-type quark multiplets 
and the down-type quark multiplets respectively and the gauge interactions do not 
mix cj) 1 and (j> 2 with the other Higgses. Thus, we can use the formulae in fl4|| , 
which are expressed by the step functions just like (jl^). Here we will use the RGEs 
that is gained by setting all 9's to one above the SUSY-breaking scale ms and all 
6's except 6q, 6q and 6^ to one below the scale ms at the RGE formulae in Q , 
Of course the correct RGEs below ms depend on the mass spectrum of the SUSY 
particles, but this difference is thought to be small and negligible. Further we will 
also neglect the r Yukawa coupling h" for simplicity. Neglecting h" does not lower 
the upper bound m% because the coefficient of h 2 in the RGE of h'[ is positive. 
Thus the RGEs we are using are 



f)h" 



above mg and 



«3 



16tt 2 

K 
16tt 2 



K 

16tt 2 

K 

16tt 2 



~gl - 3.92 - ^gi + + hi 2 



ll - -9? - —9? + -ho 2 + -h[ 2 

7q - -9? + —9? + -/l9 2 + -/i'i' 2 

« 4 W2 20 tfl 2 2 2 1 



(22) 
(23) 

(24) 
(25) 



below m s . Here g 3 = g si g 2 =g and g x = ^/bjig' . 

Next we will consider the gauge couplings. By extending the RGEs in O] in 
the same way as the above Yukawa coupling case and solving the RGEs, we obtain 
above the scale ms, 



gl 



g\ 



1 



£zMQ/M 3 y 

1 

-^ln(Q/M 2 )' 
1 



(26) 
(27) 
(28) 



and below ms, 



5 2 



<i 



gl 



constant 
1 



(for n ^ 2) 
(for n = 2) 



where Mj and M- (i — 1, 2, 3) are decided by the initial conditions: 
a 3 {m z ) =0.12, a 2 (m z ) = 0.034, a x (m z ) = 0.017, 
and the boundary conditions at Q — ms- 



(29) 

(30) 
(31) 

(32) 
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Using h% max and h'( max can be computed. 

As the criterion that the theory remains perturbative up to A, we adopt the 
condition whether one of the dimensionless coupling constants (gauge or Yukawa) 
saturates the scale A. A particular coupling A is said to "saturate" a scale A if 

^2!) < 1 (33) 
4-7T ~ 

for Q < A, and the equality in (|33| ) holds for Q = A. One can see from the RGEs 
(26)-([u]) and the conditions (|32|) that g2 is the first gauge coupling to saturate 



for every scale A. Thus one cannot set A larger than the scale A n determined by 
the saturation of gi for a given n. However, this constraint does not restrict the 
allowed region of A so strictly. For example, A 2 = 10 42 2 GeV, A 5 = 10 180 GeV, 



A 6 = 10 15 - 4 GeV and Ai = 10 10 GeV. 



5 Self-energy contribution 

The upper bound m\ discussed so far is d 2 V^' /dcj) 2 . Since the effective potential V 
is a sum of one-particle-irreducible (1PI) Feynman diagrams with zero external 
momenta, m\ is written as 

"4 = ™L- oo + n(0), (34) 

where II(p 2 ) is the self-energy of the Higgs boson. 

On the other hand, physical mass Mh is defined as a pole of the propagator, 
that is, 

M 2 =m 2 trco + n(Af 2 ). (35) 
Then from (Q) and (j3f]), the relation between ran and Mh is 

M 2 = m 2 +n(A4 2 )-n(0). (36) 

More precisely, 

Ml = m 2 h + Ren(Af^) - ReII(0). (37) 

The imaginary part of n(Af ? 2 ) — 11(0) contributes to the decay width of the Higgs 
boson, and if we suppose Mh < 2M\y, where Mw is the physical mass of the W 
boson, the Higgs is stable at tree level and (|36| ) is correct. 

This correction An(M 2 ) = n(M 2 ) — n(0) mainly comes from one-loop elec- 
troweak interaction effects. 

This correction is larger than the MSSM case due to the existence of the extra 
Yukawa couplings, but still small and lift up the upper bound mh by 2 GeV at 
most. 



6 Results 

The results including the self-energy contribution is plotted in Fig{l] in the case 
of n = 2. Since we are interested in the absolute upper bound on the mass of the 
lightest Higgs, only large tan/3 case (tan/3 = 40) is plotted. The lines are the case of 
A = 10 16 , 10 12 , 10 8 , 10 4 GeV from bottom to top. The common soft SUSY-breaking 
stop mass parameter m is set to be 1 TeV in all the cases. We can see from this 
figure that setting the saturation scale A larger, the upper bound m,h becomes more 
stringent. 



Because the RGEs (E2)-(25) have n-dependence only through the terms involving 



the weak gauge couplings g and g', the n-dependence of m,h is quite small. In fact, 
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since the n-dependence of the result becomes larger as A is larger and g becomes 
no longer small at 10 15 GeV for n = 6 (Ag = 10 15 4 GeV), we can see the most 
enhanced n-dependence of rrih when we compare the two cases: n = 2 and n = 6 
for A = 10 15 GeV. This situation is plotted in Fig.|. 

From this figure we can see that the n-dependence of nth can be neglected. 
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Xt/ms 

Figure 1: The upper bound mh in the case of n — 2. The lines corresponds to the 
case of A = 10 16 , 10 12 , 10 8 , 10 4 GeV from bottom to top. The parameters, tan/3 
and the common soft SUSY-breaking stop mass to, are set to be 40 and 1 TeV 
respectively. 
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